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Numerical Scheme 

This numerical scheme belongs to a family of codes known as UNCLE (UNsteady Com- 
putation of fieLd Equations) as reported by Whitfield (1995), that is being used to solve 
problems in a variety of areas including compressible and incompressible flows. This 
derivation is specifically developed for general unsteady multi-blade— row turboma- 
chinery problems. 

The scheme solves the Reynolds-averaged N-S equations with the Baldwin-Lomax 
turbulence model. In Chen’s work (1993), the governing equations are cast in time-de- 
pendent curvilinear coordinates with conservative variables written in inertial Carte- 
sian coordinates, 

< dF , dG , dH _ dG d , dH d m 

dr d£ drj drj d£ 

Q is the conservative variable vector, F,G,H are the convective flux vectors in the curvi- 
linear directions £/ V> £ respectively. G<l and H d are the viscous flux vectors. A finite vol- 
ume discretization is applied to the above, and a general implicit scheme is used to inte- 
grate the discretized equation in time with second order accuracy as explained in 
Janus (1990). 

The resulting numerical expression can be symbolically written as a nonlinear system 

S ( Q" +1 ) = 0 (2) 

where n+1 denotes the new time step. The convective flux in Eq. (2) is evaluated by 
Roe’s flux-difference-splitting which is first order accurate in space. High order accu- 
racy of the flux is obtained through the use of limiters. The viscous flux is central-dif- 
ferenced and is second order accurate. 
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For a flow which is steady in time Q" +1 will converge to a value which is independent 
of n. If one is not interested in the transient state, all that matters is the value of Q" +1 
as n becomes larger. However for a time dependent flow, the intermediate values of Q n+1 
are of importance for they represent the state of the fluid at time n+1. For unsteady 
flows, Equation (2) must be solved to within truncation error at each instant in time. 
The UNCLE code uses the Newton procedure, 

S'(Q n+u ) (Q n+U+1 - Q n+U ) = S ( Q n+1Jc ) (3) 

to solve Eq. (2). The variable k is the Newton iteration count and $'( Q n+U ) is the Jaco- 
bian matrix. The inviscid contributions to S' are based on a flux-vector splitting tech- 
nique. Following the work of Taylor (1992), the viscous contributions are evaluated 
numerically. The inversion of Eq. (3) at any k is by means of a symmetrical Gauss-Sei- 
del iteration procedure. 

Boundary Conditions 

The inlet boundary is specified by a characteristic-based one— dimensional boundary 
condition in an attempt to preserve the specified radial distribution of total tempera- 
ture and total pressure. The exit boundary is also specified by a characteristic type 
boundary condition making use of the radial momentum equilibrium that accounts for 
swirl produced downstream of the rotating blades. Adiabatic-wall and no-slip condi- 
tions are employed at impermeable surfaces. For boundaries at the interface of adja- 
cent blade passages, the periodic boundary condition is used. 

Tip Treatment 

Tip clearance is modeled by assuming the gap region is periodic. This means the flow 
is transported tangentially without a loss of mass, momentum and energy. The effect 
of vena contracta is considered by reducing the size of the real tip clearance to half. De- 
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tail of this model can be found in Kirtley at al. (1990). The present tip treatment is not 
exact because the tip model leads to a simplified grid that has a void inside the tip re- 
gion. It is the authors’ feeling that the correct tip treatment could be important to re- 
solve the wake close to the casing. An effort of computing instead of modeling the tip 
region by adding another grid block inside the tip region is underway. 

Thrbulence Model 

The turbulence model used in this work is the Baldwin-Lomax model. A wall function 
approximation has been incorporated into the code which allows the center of the first 
grid cell off the wall to be positioned much further out in the flow, hence reducing the 
total number of grid points. 

Results and Discussions 

An H-grid of 132 points stream wise, 51 points spanwise and 41 points pitchwise was 
used. The original computed pressure ratio vs. mass flow is shown in Fig. 1. These nu- 
merical results were obtained by time averaging the time accurate simulations, the 
time step being l/40th of the blade passing cycle. As a result, each operating point of 
Fig.l requires 10.4 CPU hours on Cray C90. 

With the exception of the pressure ratio and efficiency vs. mass flow plots, all of the 
remaining results that were presented at the ASME blind test case session are for an 
arbitrary instant in time. Subsequently, it was found that the results were time depen- 
dent. It appears in further examination of the flow field at the maximum flow condition 
that both the near hub and tip regions are the sources of the unsteady behavior. The 
time average spanwise total pressure ratio distribution at station 4 (postdiction) is 
shown in Fig. 2 along with the predicted results at an arbitrary instant in time. The 
time average results appear to be in better agreement with the measurements near the 
hub region than the predicted results. The postdicted result of Fig. 2 is the time average 
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over a cycle of rotor revolution. Further study is needed to assess the cause of the un- 
steady behavior. 
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Figure 1 Pressure Ratio vs. Mass Flow 
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Figure 2 Spanwise Distribution of Total Pressure Ratio at Station 4, High Flow. 
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